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, A simple proof is given for a generalized form of a theorem of Soshnikov. The 

£NJ ' latter states that the Janossy densities for multilevel determinantal ensembles sup- 

ported on measurable subspaces of a set of measure spaces are constructed by du- 
alization of bases on dual pairs of ^-dimensional function spaces with respect to a 
pairing given by integration on the complements of the given measurable subspaces. 
The generalization extends this to dualization with respect to measures modified by 
| arbitrary sets of weight functions. 

m 
o 

^ 1. Multilevel determinantal ensembles. 

Let {(Tj, d/j,j)}j = i... m be a set of measure spaces and {Hj := L (Tj, dfj,j)}j = i... m the 
Hilbert spaces of square integrable functions on them. Let {/ a }a=i..JV an d {h a } a= i,,,N be 
bases for a pair of iV-dimensional subspaces Hi C Hi, H m C H m , respectively. Suppose we 
are given m — 1 functions {gj+x,j}j=x...m-x on the product spaces Ij+i x Tj, such that the 
corresponding integral operators 



9j+x,j ■ Hj — > H j+ i 



9j+x,j{f){xj+x) ■= / Qj+iAxj+ux^fix^dfijixj), 

together with their transposes 

9j+x,j '■ Hj+i — > H j 
9j+i,j{f)(xj) := / g j+1 , j (x j+ i,x j )f(x j+ i)dfi j (x j+1 ), 



(1.1) 



(1.2) 



are well defined injective maps on a sequence of dense subspaces Hi C H i , H% = g2x{Hi) C 
H 2 , . . . , H m = g m , m -x{H m _i) C H m , as are their composites: 

9kj ■= gk,k-x ° . . . o gj+ij. (1.3) 
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We assume that Hi C -Hi, H m C .H m , and denote the respective images as -H2 := #21 (Hi), . . . , 
#m := 9m,m-i{Hm-i) and ff"*" 1 := <4, m _i(# m ), . . . , H 1 := g^H 2 ). Assuming further- 
more that the N x N matrix: 

A ab := f a (x 1 )g* rnl {h h ){x 1 )d^ 1 {x 1 )=j m i(/ o )(x m )(/i&)(x m )d/i m (x m ) (1.4) 

is nonsingular, it follows that the pairs of spaces {Hj,H J } may be viewed as mutually dual, 
and identified through the d\x 3 - integration pairings. Making a PLU decomposition of A 

A := PLU, (1.5) 

where L and £7 are, respectively, lower and upper triangular matrices, normalized, e.g., with 
equal diagonal entries (unique up to 2 N sign ambiguities on the diagonal), and P a permutation 
matrix, we may form bases: 

N N 

} == E( PL )a6/* 0i m) == E ^ (1-6) 

6=1 6=1 

for the spaces Hi and .H m , respectively, as well as for the sequence of dual spaces {Hj,H j } 
through composition with g 3 \ and g^: 

{V# } := 9iiM } )}a=i...iv, {<tf } := ^(0i m) )}a=i...iv. (1.7) 
These are, by construction, mutually dual, 

/ ^\x 3 )^\x 3 )d H {x 3 ) = 5 ab . (1.8) 

An example of the above construction consists of choosing the sets V j to be intervals on 
the real line, with Lebesgue measure, and the functions f a , h a and {gjj-i} of the form: 

f a {xi) := xr X e-* v ^\ h a (x m ) := x^e~^ x ™\ 

g ■_ 1 (x-,x--i) := e Xi - lXi ~^ Vi -^ Xi -^ +v ^ x ^\ 

where the V 3 (xj)'s are suitably defined functions for which the integrals involved are convergent 
(e.g., polynomials of even degree with real, positive leading coefficients). This case arises in 
the study of random multimatrix chains models [EM, BEH1, BEH2]. The reduced probability 
density for eigenvalues {xa}j=i... m ,a=i...N is given by the determinantal formula 



m — l 



Pn(^) = ^-det(tf)(x< 1) ))det(^)(x^ ) )) ]J detfo i+1J (*? +1) > x[ j) )) (1.10) 
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where the partition function Zjv,m is defined so as to normalize this to a probability measure. 

More generally, let us suppose that the functions f a , h a , gjj-i are chosen so that the 
expression (1.10) defines a probablity measure on the ensemble YlT=i(^j) N ■ Such ensembles 
are sometimes referred to as "determinantal ensembles" [BS, S]; other examples include, e.g., 
"polynuclear growth" models [PS, J]. 

Define the following set of functions on the product spaces T; x : 

N 

K ij (z i ,z j ):=Y l 1>!? ) (xi)<l>PM 
0=1 

and 

■Kij(Xi, Xj) . Kij^Xi, Xj) 9ij(Xii Xj) , (1.12) 

where gij(xi,Xj) := if % < j. It may then be shown [EM] that that the probability density 
(1.10) may equivalently be expressed in the form: 

PN, m (x^) = det^sW,*^)) , (1.13) 

where Kij(xa \ x£ ) is viewed as the ((i, a), (j, b)) element of a matrix of dimension Nm x iVm, 
labelled by pairs of double indices 1 < i,j < m, 1 < a, b < N. By integrating (1.13) over 
a part of the variables, it follows [EM] that the correlation function giving the probability 
density for finding kj elements in Tj at the points {xf , . . -x^} for j = 1 . . .m is similary 
expressed by the X^jli x Sjli determinant 

Pfc!,..*™^^...*^^ • (1-14) 

The functions K^- may also be viewed as kernels of integral operators 



Kij : Hj — > Hi 



Kij(f)(xi) := / K lJ (x l ,Xj)f(x j )dfi j (x J ) 



(1.15) 



which, again, map the finite dimensional spaces {Hj} to each other. Note that the various 
ify's may all be obtained from K\ m by composition on the left and right with operators gif 

= gn o K lm o g mj , (1.16) 

and that, in particular, adjacent ones are related by 

Kij = g^i-i o Ki- ly j = K iJ+1 o g j+ i yj . (1-17) 
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It follows from (1.8) that, when restricted to Hj and H\ Kjj and K*- act as identity operators 

K JJ ^) = ^, ipeHj, K*^) = <f>, 4>eW. (1.18) 

In the following, we denote by K, g and K the m x m matrices of integral operators 
acting on the direct sum 

H:=®f =1 Hj, (1.19) 

with matrix entries K^, gij and Kij acting on the component spaces Hj. 

Let {Ij C Tj}j = i mm _ m be measurable subsets of the spaces Tj, and let Xj denote the 
characteristic function of Ij C Tj. Denote by xi the direct sum of these functions, viewed as 
an operator on H, acting by multiplication by the various x/s on the component spaces Hj's, 
and assume that this leaves H invariant. Let 

K Xl :=Ko X i (1.20) 

denote the composition of these operators, also a matrix integral operator, and let 

R Xl := (1 — K Xl )~ 1 o K Xl (1.21) 

be its Fredholm resolvent. Then the matrix components of this resolvent, denoted Rij, are 
operators with integral kernels Rij(xi,Xj). These also define certain correlation functions, on 
the product YiT=i if ■> the so-called Janossy densities, given by a formula similar to (1.14), 
namely: 

PiL.fcJ^V--^^ , (1-22) 

where the normalization constant Cf' m is defined to be the Fredholm determinant: 

Cf' m :=det(I-K Xl ) (1.23) 



(which equals the probablity, under the original distribution (1.10), of having no elements 

i=i L i 



within the subset / = fl^Li if)- The correlation functions (1.22) give the probability density 



of finding, for j = 1 . . .m, exactly kj elements in Ij at the points {xf \ . . . x k f }j=i... m • 

A theorem of Soshnikov [S] expresses this distribution in terms of a new set of functions 
{i/jti') , <f>^\gij}, analogous to 4> i9ij\- To define these functions, one begins by defining, 

as in (1.3), integral operators 

9kj ■= 9k,k-i °!c ■■■ Ojc 9j+i,j (1-24) 

k — 1 j + 1 
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where the symbol o jC denotes composition of integral operators by integration only over the 

3 

domain Ij complementary to Ij within Fj. (Note that gj+i,j = Qj+i.j-) Analogously to (1.4), 
we define the pairing matrix 

A ib '■= / fa{xi)g* rnl {h h ){x 1 )dn 1 (x 1 ) = / g m i(fa)(x m )(hb)(x m )dfi m (x m ), (1.25) 

1 m 

and again assume it to be nonsingular. Once again, making a PLU decomposition of A 1 

A 1 := P 7 Z/[/ 7 , (1.26) 

we may form new bases 

N N 

> := E^UVfc ^ == Y,^aha (1-27) 

6=1 a=l 

for the spaces Hi and H m , as well as bases for the sequence of dual pairs of spaces obtained 
by composition with gji and g^j'- 

A j) ■= 9ji(A 1] ), S' } ■= T m 0r ] )- (1-28) 
These are mutually dual in the sense that 

^\x j )$\x j )d H {x j ) = 5 ab . (1.29) 



Define, as before, a set of integral kernels: 

N 



KyixuXj) := 'EWMtPM (1-30) 

a=l 

and 

Kij^Xi, Xj^) := Kij^Xi, Xj) gij(xi,Xj), (1.31) 

where again, gij(xi,Xj) := if % < j, and denote by Kij, and the corresponding 
integral operators mapping Hj — > Hi. The associated m x m matrix integral operators acting 
on H are similarly denoted K, g and K, with 

K:=K-g. (1.32) 

Then Soshnikov's theorem states that the resolvent operator R XI entering in the definition of 
the Janossy distribution coincides with the operator K XI := K o xi so constructed. 
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Theorem 1: 

R XI =K XI . (1.33) 

A proof of this theorem is given in [S] , based on decomposition of the spaces into various 
subspaces, and an inductive construction, but it is rather long and intricate. In the following 
section, a very simple direct proof will be given. 

Since no further effort is required, we will actually prove an obvious generalization that 
reduces to the above result as a special case. Namely, instead of choosing measurable subsets 
Ij C Tj and their characteristic functions Xji we ma y replace the latter by any set of square 
integrable functions {wj G L 2 (Tj, dfj,j)}j=i... m , and generalize the definition of the matrix A 1 
and the functions {^\ (f>^\gij} accordingly to: 



Kb '■= / fa{xi)g* m i{hb){xi){l - wi(zi))<W:ri) 

= / 5ml(/«)M(^)(Sm)(l - W m (x m ))d/i m (l m ), 

Jr m 

with PLU decomposition 



(1.34) 



A w ._ p w L u, UW ^ ( L35 ) 

where, denoting by o v . the operation of mutiplication by a function Vj in L 2 (Tj, d^ij) followed 
by composition, we replace the definition of gk,j in (1.24) by 

9kj ■= 9k,k-i oi-itfk-! o . . . o 1 _ Wj+1 og j+1J . (1.36) 

Then 

N N 

A 1] ■= Y.( PWL ^abfa ^ ■= Ts^baha, (1-37) 

a=l b=l 

define bases for the spaces Hi and H m , respectively, while bases for the sequence of dual 
spaces obtained by composition with Tjj i and k , as defined in (1.36), are again given by 
(1.28), with the dualization pairing given by 



^(x^ix^l - WjixiMniixj) = S ab . (1.38) 



Replacing the characteristic functions xi = (Xi ■ ■ ■ Xm) by the set of weight functions w := 
(wi, . . . w m ), the definitions of the integral operators and kernels K, g, K retain the same form 
and the theorem becomes 

Theorem 2: 

R w = K w := K ow. (1.39) 
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where 

R w := (1 — K w )~ l o K w (1.40) 

is the Fredholm resolvent of 

K w :=Kow. (1.41) 

Such a generalization is of use in multimatrix models, since it allows us to replace, for 
example, the characteristic functions Xj on subintervals Ij C Tj by weighted characteristic 
functions Yla=i K a,jXi a ■ over unions of disjoint subintervals Ij = U^ J =1 / Qi j. The correspond- 
ing Fredholm determinant (1.23), expanded as a power series in the coefficients K a j becomes 
a generating function for the probabilities of having given numbers of eigenvalues within the 
subintervals, as in the 1-matrix case [TW]. 

2. Proof of the theorem. 

The equality (1.39) may equivalently be written as 

K w oK w = K w -K w , (2.1) 

or, explicitly, 

m 

^ Kij o w . K jk = K ik - K ik (2.2) 

3=1 

(where composition on the right by w is omitted, since the equality will be shown to hold 
without it). Substituting (1.12) and (1.32) in (2.1), this is equivalent to 

Ko w k-go w K-Ko w g + go w g + g- g = K- K. (2.3) 

This relation follows as a consequence of four identities relating the various summands on the 
left, each of which is easily proved: 

g° w g + g-g = o (2.4) 
(K o w K)ij = g a o K Xj + 5nK l3 - K im o 1 _ Wm g mj - K irn S mj (2.5) 
(g o w K)ij = ga o K Xj + 5nKi 3 - Kij (2.6) 

\K O w g^ij -Kij Kim °1 — w m gmj Ki m 5 m j (2.7) 

Proof of (2.4): From (1.3) and (1.36) we have, for k - 1 > j > % + 1, 

gkj ° W j gji = gkj ° gji - gkj °i-wj g 3 i 

= gkj ° gji — gkj+i ° gj+i,j °i—w j gji (2.8) 
= gkj gji - 9k,j+i 9j+i,i, 



(2.9) 



(2.10) 
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while for j = k — 1, 

9k,k-i °^ fe _i 9k-i,i = 9k,k-i ° 9k-i,i — 9k,k-i °i-w k -i 9k-i,i 
= 9k,k-i ° 9k-i,i - 9ku 

and for j = i + 1 

g/c.i+l <7i+l,i = ° 9i+l,i — 9k,i+l 9i+l,i 

= 9ki - 9k,i+2 ° 9i+2,i- 
Summing over j and cancelling all intermediate terms gives the result 

k-i 

{g ° w g)ki = 9k i ° w i to* = 9ki ~ 9ki - i 2 - 11 ) 
j=i+i 

(Note that this means that the integral operator g w is the Fredholm resolvent of g w .) 
Proof of (2.5): For j ^ m, we have 

Kij o w . Kjk = Kij o Kjk — Kij °i- W j Kjk 

= o Kjk - Ki,j+i ° 9j+i,j °i-wj Kjk (2.12) 
= K^ o Kjk - K iJ+1 o K j+ljk . 

Therefore, summing and cancelling the intermediate terms gives 

m 

^ o w . Kj k = Kn o Kik - K irn °i-wj K mk . (2-13) 
j'=i 

By (1.8) and (1.38), the operators K\\ and (K mm o (1 — w m ))* act as the identity on the 
spaces Hi and H m , respectively, and hence 

Kn o K lk = K lk , K im o 1 _ Wm K mm = K im . (2.14) 

It follows from (1.16) and the corresponding relations 

Kij = 9n oi_ Wl Ki m °l-w m 9mj (2-15) 

for the ify's that, for i ^ 1, j ^ m, 

K a O Kik = 9il ° K ik , K im °l-w m K m j = Kim °l-w m 9mj ■ ( 2 -16) 

Combining these relations with (2.13) leads to (2.5). 
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(2.18) 



Proof of (2.6): For j < i — 1, we have 

9ij °wj Kj k = gij o K jk - g t j oi- Wj K jk 

= o K jk - g iJ+1 o g j+1J o x _ w . K jk (2.17) 
= g^ o Kjk - 9i,j+i ° Kj+i,k, 

while for j = i — 1, 

S'i.i-i u,„i Ki-i,k = 9i,i-i ° ^i-\,k — 9i,i-i i-wi_i ^i-\,k 

= 9i,i-i ° ^i-i,fc - -fQfc, 
Again, summing, and cancelling the intermediate terms gives: 

m j — 1 

o Wj . iT jfc = g-y o Wj . K ifc = g n o Ki fc - K ik , (2.19) 
i=i i=i 
which is (2.6) (the case i = 1 being trivially satisfied). 

Proof of (2.7): For m > j > k + 1, we have 

Kij o w . g jk = Kij o g jk - Kij o t _ w . g jk 

= K^ o g jk - K iJ+1 o g j+ ij °i- Wj 9jk (2.20) 

= Kij ° 9jk - ° 9j+i,k, 

while for j = k + 1, 

Ki,k+i °w k+1 9k+i,k = Ki,k+i ° 9k+i,k — Ki,k+i °i-w k+1 9k+i,k 
= K ik — K ijk+2 o g k+ 2,k, 
Again, summing and cancelling the intermediate terms gives: 



(2.21) 



^ Kij o w . g jk = ^ K^ o w . g jk = K ik - K im o x _ Wm g mk , (2.22) 

j=l j=k+l 

which is (2.7) (the case j = m again being trivially satisfied). 

Combining the four relations (2.4), (2.5), (2.6) and (2.7) gives (2.3), hence proving the 
theorem. 

Remark. Although this theorem has been formulated with respect to probability measures 
and Hilbert spaces of functions on them, it may clearly be extended to more general cases 
involving, e.g., complex measures, provided all the integrals appearing are convergent, and 
the related integral operators and their composites are well-defined. 
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